Helical rotating turbulence. Part II. Intermittency, scale invariance and structures 
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We study the intermittency properties of the energy and hehcity cascades in two 1536'^ direct 
numerical simulations of helical rotating turbulence. Symmetric and anti-symmetric velocity incre- 
ments are examined, as well as probability density functions of the velocity field and of the helicity 
density. It is found that the direct cascade of energy to small scales is scale invariant and non- 
intermittent, whereas the direct cascade of helicity is highly intermittent. Furthermore, the study 
of structure functions of different orders allows us to identify a recovery of isotropy of strong events 
at very small scales in the flow. Finally, we observe the juxtaposition in space of strong laminar and 
persistent helical columns next to time-varying vortex tangles, the former being associated with the 
self-similarity of energy and the latter with the intermittency of helicity. 



I. INTRODUCTION 

Turbulence is often referred to as the last unresolved 
main problem of classical physics. The diversity of ap- 
plications of turbulent flows, from geophysics and astro- 
physics to engineering, and the observed complexity and 
lack of predictability of such flows both make the problem 
difficult to tackle, or even to define. Part of this complex- 
ity is due to the fact that turbulence comes in intermit- 
tent "gusts," strong events that are scarce, but still more 
frequent that what could be expected if normally dis- 
tributed. These gusts give rise to the well-known break 
down of scale invariance in the flow 

Intermittency is a highly spatially and temporally lo- 
calized phenomenon. It is believed to be associated only 
with a forward cascade of an ideal invariant (a quan- 
tity conserved in the inviscid case), corresponding to the 
transfer of this quantity towards small scales with con- 
stant fiux as a result of the non-linear coupling between 
modes; it is not directly related with the dimensional- 
ity of the problem. As an example, the energy cascade 
from larger to smaller scales in isotropic and homoge- 
neous three dimensional turbulence is intermittent 
Magnetofluids, i.e., conducting fluids where the velocity 
field is coupled to the magnetic field, are intermittent in 
two 3 and in three dimensions and in both cases a 
direct cascade of energy takes place. These examples are 
in contrast with the two dimensional neutral fluid case, 
for which the conservation of vorticity leads to an inverse 
energy cascade (a transfer of energy to the large scales 
with constant flux). While the direct cascade of enstro- 
phy in this case is intermittent, there is evidence that the 
inverse cascade of energy is scale invariant and probabil- 
ity density functions of velocity fluctuations are close to 
Gaussian For the case of rotating turbulence, stud- 
ies of intermittency are not numerous since they require 
high Reyolds numbers together with low Rossby num- 
bers; moreover, since in this case energy undergoes both 
a direct and an inverse cascade @, 0j @1 1 unclear a pri- 



ori whether intermittency is to be expected in the small 
scales or not, although, until the present work (see Sect. 
IIV[) , there was no documented example of self-similar di- 
rect energy cascade to our knowledge. In experiments of 
freely decaying rotating turbulence [9j , it was found that 
rotation signiflcantly decreases intermittency, although 
strictly scale-invariant (non-intermittent) statistics were 
not found. Such experimental results are also in agree- 
ment with recent analyses of data stemming from direct 
numerical simulations (DNS) of non-helical rotating tur- 
bulence El • Other laboratory experiments re- 
ported an even more pronounced reduction of the inter- 
mittency when rotation is present. 

Although intermittency is believed to take place at 
small scales, strong events can affect the dynamics of the 
large scales, specially in systems close to criticality. As 
an example, it was shown in [l^, [ij, [l^ that local fluc- 
tuations of the kinetic helicity Hy = u • V x u with u the 
velocity, can explain phase and amplitude variations of 
the 22-years solar cycle. Also, intermittency is known to 
affect the transport of momentum in atmospheric surface 
layers [l6j . 

Considering the large number of degrees of freedom 
typical in turbulence (typical mesoscale flows in the at- 
mosphere of the Earth can have, for instance, about 10^^ 
degrees of freedom), many of the approaches to turbu- 
lence rely on an assumed scale invariance to model the 
flow statistics at small (often unresolved) scales. The 
search for self-similar quantities in three-dimensional tur- 
bulence is a long-standing problem, which would re- 
late its study with critical phenomena and the out-of- 
equilibrium statistics of systems with a large number of 
modes, and which would allow the use of tools from quan- 
tum field theory, condensed matter, and statistical me- 
chanics. 

In a previous paper, hereafter referred to as Paper I 
we presented results from two massive numerical 
simulations of helical rotating turbulence. Rotation and 
helicity arc present in many atmospheric phenomena, as 
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e.g., supcrccU storms and tornadoes [l8|, ll9|, |20[. The 
simulations confirmed the scaUng laws for the energy and 
helicity spectra predicted in [2l|. The development of 
anisotropics in the flow, as well as scaling laws in the di- 
rections parallel and perpendicular to the axis of rotation 
were also studied, considering both global measurements 
of anisotropy as well as spectral quantities. In this paper 
we study the intermittent properties of both the direct 
cascade of energy and the direct cascade of helicity in 
helical rotating turbulence, using the data from the sim- 
ulations presented in Paper I. A decomposition into di- 
rections parallel and perpendicular to the axis of rotation 
is used to study anisotropy. Details of this decomposi- 
tion, as well as the definitions of increments, structure 
functions, and scaling exponents used to quantify the in- 
tcrmittcncy are given in Sect. |TT1 The results of the de- 
composition performed on the data stemming from the 
high resolution DNS are presented in Sect. IIIIl where the 
possible recovery of isotropy at small scales is discussed. 
Averaging over different directions and over time, up to 
4.3 X 10"'^^ data points were used to obtain statistical evi- 
dence that the direct cascade of energy is scale invariant, 
while the direct cascade of helicity is not, as discussed in 
Sect. IIVI This result is obtained from the study of struc- 
ture functions as well as probability density functions of 
velocity and helicity increments (Sect. |V|. Finally, the 
development of structures at large and small scales in 
the flow is discussed in Sect. IVII and the conclusions are 
presented in Sect. IVIII 



As a result, if the direct cascade is scale invariant, for an 
isotropic and homogeneous flow Sjj{£) (p^^^, and depar- 
tures of the exponents from this relation are a signature 
of intermittency. 

The study of intermittency in the direct cascade of he- 
licity has been less explored. In the isotropic and homo- 
geneous case, one can study it using for example struc- 
ture functions based on the helicity flux, which is a third 
order field [2^, [2^. In real space, this flux can be written 
in two different ways, which follow from the r.h.s. of the 
Karman-Howarth theorem for the helicity [2^, [2^ : 



u(x) 



£ 



(u(x) X u(x + £)) 
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Sf, (3) 



or in terms of structure functions of the velocity and vor- 
ticity nil. 



-i(fc(x,^) (<5u(x,^))' 



-I- 



(4) 



where e is the helicity injection rate. Both relations are 
equivalent, as it is easy to see from Suj ^ Su/i How- 
ever, here we want structure functions that can be associ- 
ated to the spectral scaling, and must therefore be based 
on second order quantities. There are two candidates 
that follow from the expressions in the time derivative 
appearing in the Karman-Howarth theorem for helicity 
|2i,|2i: 



II. VELOCITY AND HELICITY INCREMENTS 



([u,(x + £)- u,{^)] [u,i^ + £)- u,{^)]) , (5) 



A. Increments 

To study intermittency in the direct cascade of energy, 
we will use longitudinal increments of the velocity field 



Mx,^) = [u(x + £)-u(x)]. Y, (1) 

where the increment £ can be in any direction. Structure 
functions of order p are then defined as 



Spi£)^{Su''i^,£)), 



(2) 



where the brackets denote spatial average over all values 
of x. The structure functions depend on the direction 
of the increment, and no assumption about isotropy or 
axisymmetry has yet been made. 

The Sp structure functions arc of interest because for 
p = 2 they are related to the two-point correlation func- 
tion of the velocity, and thus to the energy spectrum 
through the convolution theorem. Therefore, a power 
law behavior can be expected for scales corresponding to 
the inertial range. Also, for the isotropic and homoge- 
neous case, 5*3 is related to the energy flux and scales 
linearly with the increment £ in the inertial range [2^ . 



(5u(x,£)-(5w(x,£)) 



(6) 



Both quantities are Galilean invariant, as is required for 
the structure functions to be well behaved [2^ . The first 
quantity is a second order tensor with indices i, j, while 
the second is a scalar. For simplicity, here we use the 
second expression, although both are related to the an- 
tisymmetric part of the two-point correlation tensor for 
the velocity. Then, structure functions for the helicity 
can be defined as 



np{£) = {[Su{x,£)-Su{^,£)f). 



(7) 



With this definition, for isotropic and homogeneous 
turbulence the assumption of scale invariance leads to 
7ip{£) ~ £P^'^. Note that Hi is second-order in the veloc- 
ity whereas Si is first order. 



B. Parallel and perpendicular directions 

The development of anisotropics in a rotating flow has 
been studied in experiments [29l . [sol . [3ll | and in numerical 
simulations (see e.g, [s^V Anisotropy was shown to de- 
velop in a range of Rossby numbers such that non-linear 
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FIG. 1; The 12 generators used to compute increments in 
the x-y plane, and the generator in the z direction. The 
crossings ol dotted lines indicate grid points in the numerical 
simulation. 



interactions are not completely damped with the scram- 
bling effect of inertial waves [3^ . The presence of rota- 
tion breaks down the isotropy of the flow, introducing a 
preferred direction. Energy, as a result of resonant triad 
interactions, is transferred preferentially towards modes 
in spectral space perpendicular to the axis of rotation 
@j 0, HI] 1 although linear effects may also be relevant in 
the formation of the structures 0|. 

In isotropic and homogeneous turbulence, it is a com- 
mon practice to study velocity increments (or two-point 
correlation functions) in only a few directions, or to av- 
erage over different directions as in, e.g., the leading 
(isotropic) term of the SO (3) decomposition [H, Is^ . 
Given the preferred direction in our problem, and the 
natural axisymmetry associated with it, we will be in- 
terested in increments parallel and perpendicular to the 
angular rotation $7, which will be denoted respectively 
as i\\ and In principle, in the plane perpendicular to 
rj, the increments can be taken in every possible di- 
rection and later averaged to obtain structure functions 
that only depend on the scalar increment l^. However, 
this requires interpolation of the three components of the 
velocity in the three dimensional space every time an in- 
crement docs not reside on a computational grid point. 
To avoid the computational cost of this interpolation, we 
follow the procedure introduced in [IJ for the isotropic 
case. We only compute increments given by the product 
of an integer times a generator vector, with the generator 
vectors chosen such that they fill as uniformly as possi- 
ble circles in the plane perpendicular to fJ, and such that 



their product with an integer always falls on a grid point. 

Twelve generators were used in the x-y plane: (1, 0, 0), 
(1,1,0), (2,1,0), (3,1,0), (0,1,0), (-1,1,0), (1,2,0), 
(-2,1,0), (-1,2, 0), (1,3,0), (-3,1,0), and (-1,3,0) (in 
units of grid points in the simulation). These genera- 
tors, plus the 12 generators obtained by multiplying them 
by —1 (or equivalently, considering negative increments) 
cover the plane in an approximately uniform way (see 
Fig. [T]). In the z direction (the direction parallel to 
the axis of rotation) the generator for the increments is 
the vector (0,0,1). Given these generators, increments 
in Eqs. ^ and ([7]) are created by multiplying the gen- 
erators by integer numbers. With these choices, all in- 
crements reside on grid points and no interpolation in 
the computation of the structure functions for each di- 
rection is required. Once structure functions for all di- 
rections have been computed, structure functions in the 
perpendicular direction Sp(l\_) and TLp{l±_) are obtained 
by averaging over the results for the 12 directions in the 
x-y plane. As the generators have different lengths, in- 
terpolation in this step is required, but it is less costly 
from the computational point of view as only interpola- 
tion of scalar one-dimensional functions is needed. The 
structure functions in the parallel direction, Sp{l\\) and 
Tlp{l\\), are obtained directly from the generator in the 
z direction. Finally, average in time (using snapshots 
of the velocity field at different turnover times) can be 
computed. 

In our case, we use for the run with the larger rotation 
rate (0 = 9), ten snapshots of the velocity field spanning 
ten turnover times, from i = 20 to 30 (see paper I). As a 
result, considering that each snapshot has 1536'^ « 3.6 x 
10^ grid points, and considering the twelve generators 
used, each increment in Sp{l±) and Ti.p{l±) results from 
an averaging over 4.3 x 10^^ data points. In the case 
of 5^(^11) and Hp{l\\), each increment is obtained using 
w 3.6 X 10^° data points. 

III. STRUCTURE FUNCTIONS 
A. Velocity structure functions 

Two simulations were used for the analysis, described 
in more detail in Paper I. One of the simulations (here- 
after, run A), has il = 0.06, a Rossby number Ro w 8.5, 
and was continued for almost 10 turnover times. The 
other simulation (run B) has = 9, a Rossby number 
Ro « 0.06, and was continued for 30 turnover times. 
Both runs have a Reynolds number Re « 5100. 

Figure [2] shows the result of computing the velocity 
structure functions in all directions for one snapshot of 
the field (at t = 30 in run B), and of averaging over 
the different directions to obtain 5*2 (^±) and S2{i\\)- At 
small scales, both structure functions scale as P, as can 
be expected for a well-resolved smooth flow in the dis- 
sipative range. At intermediate scales, an inertial range 
with power law-scaling can be identified in S2{i±), but 
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FIG. 2: Second order structure functions 52 at f = 30 in 
run B with Ro = 0.06. The dotted lines indicate the dif- 
ferent structure functions in the twelve directions given by 
the generators in the x-y plane, and the thick solid curve is 
the average ^(^x). The thick dashed curve corresponds to 
increments in the z direction and is 52(^11). 
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FIG. 4: Sixth order structure functions Sq{Ia_) and S6(^||) in 
run B at different times, between t = 2Q and 30. A dissipative 
range scaling ^ £^ is indicated at small scales, and two average 
slopes are indicated in the inertial ranges (see text). Note that 
the perpendicular part of the structure function dominates the 
parallel one at all scales. 
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FIG. 3: Second order structure functions S2{i±) and 52 (^||) in 
run B at different times, between t = 20 and 30. A dissipative 
range scaling ~ is indicated at small scales, and the average 
slope (2 ~ 1.41 is indicated in the inertial range. 



not in 52(^11). Indeed, S2{i\\) is smaller than S2{i±) at all 
scales (specially so at the largest scales), and shows no 
clear scaling. This is consistent with the results obtained 
in Paper I from the energy spectrum: while the energy 
spectrum in perpendicular wave vectors shows an iner- 
tial range with power law behavior and approximately 
constant flux, the energy spectrum in the parallel direc- 
tion shows no clear scaling and its associated flux decays 
rapidly with scale. 

Figure [3] shows the structure functions 6*2 (^_l) and 



S2{i\\) at four different times between t = 20 and 30 
in run B. No substantial differences are observed at the 
different times. The average scaling exponent in the iner- 
tial range S2{i±) ^i', obtained after averaging over the 
ten turnover times, is (2 = 1-41 ± 0.02, and is indicated 
in the figure by the straight line. Errors here and in the 
following are defined as the standard mean error 



\ 



N 



(8) 



where is the slope obtained from a least square fit 
for each snapshot i, N is the number of snapshots, and 

is the mean value averaged over all snapshots. The 
error in the least square determination of the slope for 
each snapshot is much smaller than this standard mean 
error for the averaged exponent. Extended self-similarity 
is not used to obtain the slopes. 

The value obtained for ^2 is in good agreement with 
predictions for helical rotating turbulence [2l|. In a ro- 
tating flow with maximal helicity, E ~ fc~^'^, which leads 
to 5*2 (£±) ~ ^Y^- For flows with non-maximal helicity the 
(2 exponent is, according to [2l|, between 1 and 1.5, with 
the value of 1 corresponding to the non-helical case. Note 
that in numerical simulations of non-helical rotating tur- 
bulence S'2(^j_) ^ i± was reported in [lol. [Tl|. 

From the behavior of the amplitude of <5'2(^||) with 
scale, it seems isotropy could be recovered at small scales 
in a rotating flow if the inertial range is wide enough for 
S'2(^||) and S2{£±) to collapse (or, in other words, for the 
effect of rotation to be negligible at the smallest scales). 
Indeed, for small values of i the parallel and perpen- 
dicular structure functions become closer. This is more 
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FIG. 5: Second-order helicity structure functions H2{£±) and 
H2{£\\) (see Eq. (0) in run B with Ro = 0.06. at different 
times between t = 20 and 30. The dissipative range scales 
as ~ £'^, consistent with the fact that H2 is quartic in the 
velocity; the average slope is indicated for the inertial range. 



evident in velocity structure functions of higher order. 
As an example, Figure [J] shows the sixth order veloc- 
ity structure function. While an inertial range in the 
perpendicular direction is still visible (the mean slope 
averaged over ten turnover times is indicated as a ref- 
erence), and the structure functions in this direction do 
not change much between different snapshots, such is not 
the case at smaller scales, or in the parallel direction. In 
the parallel direction, larger fluctuations between differ- 
ent times are observed. Moreover, the anisotropic iner- 
tial range is shortened as, at small scales, both Se{i±) 
and Sq{£\\) collapse to a single curve. This collapse takes 
place for scales smaller than £ w 0.05, with slight fluctu- 
ations in time; it indicates that the strongest events in 
the flow (which begin to dominate structure functions as 
the order increases) tend towards isotropy at the small- 
est scales. Note that the collapse is not observed in the 
S2 structure functions (see Fig. [3]), whereas Sp{£±) and 
Sp{£i\) get closer to each other at smaller scales as the 
order p is increased. 

This partial recovery of isotropy in the higher order 
moments of the velocity field can be understood as fol- 
lows: higher values of p in Eq. ([2]) leave only the strongest 
gradients contributing to ^^(.^j,) and 5p(.^||), as the con- 
tribution of mild gradients to the structure functions goes 
to zero as p is increased. It is only these strong events 
at small scales that arc more isotropically distributed, 
as indicated by Fig. [1] (at least at the Reynolds num- 
bers considered here). However, mild gradients are still 
anisotropic at small scales, as follows from Fig. [31 Sim- 
ulations with more spatial resolution at larger Reynolds 
numbers will be required to study if this transition to- 
wards isotropy only takes place for the strong events, or 
whether for sufficiently small scales the transition takes 



place for all orders, with a change in the spectral index 
of the flow (see e.g., [H, [H, for phenomenological 
treatments of non-helical rotating turbulence that con- 
sider the possibility of such a transition). 

Similarly to the stratified case, one can introduce a 
wavenumber at which rotation and nonlinear advection 
balance (see [38l[39|): 

(with e the energy dissipation rate) ; fco can be considered 
as the largest wavenumber where rotation effects are im- 
portant; the inverse of this wavenumber, £q = 27r/fco, is 
equivalent to the Ozmidov length in stratified turbulence, 
which separates the inertial range (at smaller scales) from 
the range dominated by buoyancy (at larger scales). It is 
rather remarkable that in our simulation £a Ki 0.04, close 
to the value of ^ « 0.05 where the transition in Sp{£±) 
takes place for large values of p (see Fig. 2]). However, 
a confirmation of this would require a parametric study 
varying the value of fi, which in DNS at the spatial res- 
olution considered here is out of reach with present day 
computers. 

At this point, a discussion about units is in order. The 
distance between grid points in our runs is 27r/1536 ~ 
0.004. This is also the smallest distance for which in- 
crements can be computed in the structure functions. 
Since the simulations are dealiased using the 2/3-rule, the 
largest wavenumber resolved is fcmax = 512 which corre- 
sponds to a length £^in = 27r/fc,„ax ~ 0.01. As a result, 
the velocity field at scales between w 0.004 and « 0.01 
must be necessarily smooth, and its structure functions 
should scale as Sp ~ £p as it is indeed the case (sec e.g.. 
Fig. [3]). The dissipation scale in the simulations is just 
slightly larger than £niin, which explains why the ^ £^ 
scaling extends a little bit beyond ^min (as required for 
the simulations to be well resolved). For practical pur- 
poses, we can estimate the dissipation scale to be be- 
tween w 0.01 and 0.02 (these values are consistent with 
estimations from the energy spectrum, shown in Paper 
I). Since the collapse of the parallel and perpendicular 
structure functions occurs near £ w 0.05, one may won- 
der if this scale is well resolved, or if the collapse results 
from numerical cut-off or viscous effects. We computed 
structure functions for simulations of forced helical and 
non-helical rotating turbulence at resolutions of 512^ grid 
points (see [ill |21| |) for which fcn > fcmax, and in that 
case no bump, collapse, or changes in the behavior of the 
structure functions at the smallest scales were observed. 
However, simulations at larger resolutions would be de- 
sirable to further confirm this result. 



B. Helicity structure functions 

Helicity structure functions arc computed in the same 
fashion as velocity structure functions (see Eq. ([7])); 
the functions H2{£±) and H2{£\\) for different times are 



6 



10' 



10" 



10" 
0.001 




10.000 



FIG. 6: Fourth-order helicity structure functions H4{£±) and 
iif4(^ll) in run B at different times, between t = 20 and 30. 
The average slope is indicated for the inertia! range. 
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FIG. 7: Scaling exponents (with error bars, see Table |l]) as 
a function of the order p, for the velocity (stars) and the 
helicity (pluses) in run A with Ro = 8.5, and for the velocity 
(triangles) and the helicity (diamonds) in run B with Ro = 
0.06. The dotted line corresponds to Kolmogorov scaling — 
p/3, and the dash line to = 0.71p, which represents the 
velocity exponents best. 



shown in Fig. O after averaging in the different direc- 
tions, ft should be noted that the H2 structure functions 
are effectively of fourth order in the fields, and as a re- 
sult the convergence of the statistics is not as good in 
this case as it was for the structure functions discussed 
in the previous section for the velocity. Also as a result 
of the higher-order dependence on the fields, the dissipa- 
tivc range scales as i'^. Moreover, helicity is not a positive 
definite quantity, and cancellations between regions with 
positive and negative alignments of the velocity and the 
vorticity can take place, resulting in larger fluctuations 
of the increments. Convergence of the statistics for all 
orders studied here was checked by computing the cuniu- 
lants for each moment (see e.g., [Zll. l42j ) . 



TABLE I: Order p, and scaling exponents for the velocity 
and for the helicity, with errors, for run A {Ro — 8.5) and 
run B (Ro = 0.06). 

p (p (run A) (Run A) Cp (run B) (Run B) 



0.37 ±0.01 
0.70 ± 0.03 
0.99 ±0.03 
1.25 ±0.04 
1.49 ± 0.04 
1.69 ±0.05 
1.88 ±0.05 



0.35 ±0.01 
0.59 ±0.01 
0.72 ±0.01 
0.76 ±0.02 



0.701 ± 0.007 
1.41 ±0.02 
2.14 ±0.02 
2.88 ±0.04 
3.61 ±0.07 
4.3 ±0.1 
4.9 ±0.2 



0.70 ± 0.03 
1.46 ±0.03 
2.27 ±0.07 
2.88 ±0.07 
3.1 ±0.1 



In the case of the helicity structure functions, fluctu- 
ations between different temporal snapshots are larger 
than for the velocity structure functions. This is mostly 
due to the fact that H2{£\\) displays strong fluctuations in 
time, and that the helicity becomes isotropic at smaller 
scales more rapidly with increasing p. However, fluc- 
tuations in the incrtial range of the perpendicular in- 
crements are smaller, and an inertial range can still be 
dcflned. For the second order, defining the scaling ex- 
ponents as Hp{i±) ~ i^^, we obtain ^2 = 1-46 ± 0.03. 
This value leads to a spectral scaling H(k) ^ ^Z^'^ 
(the prediction for a maximally-helical rotating flow is 
H{k) ~ fcj^'^). The scaling of the energy and of the he- 
licity obtained from the structure functions is consistent 
with the e + h = A rule for helical rotating flows, where e 
is the energy spectral index and h is the helicity spectral 
index [2l|; it is also consistent with the spectral indices 
measured in Paper I. 

Fluctuations become larger for larger values of p, and 
we were unable to compute structure functions for mo- 
ments with converged cumulants beyond p = A for run A 
and p = 5 for run B for the helicity. Specifically for run 
B, variations in the amplitude of TLp{i\\) increase with p, 
changing the scale where the structure functions become 
isotropic. As a result, for some snapshots no scaling in 
the perpendicular direction was observed, and as a rule of 



thumb snapshots for which Ti.p{i±) and 7ip{£\\) became 
of the same order at scales larger than € w 0.1 had to 
be discarded. Figure [6] shows the fourth order helicity 
structure functions for four snapshots that present a dis- 
cernible inertial range in £± (a total of eight snapshots 
was used to compute the average scaling exponents). 



IV. INTERMITTENCY IN THE DIRECT 
CASCADES 



In a self-similar flow, scaling exponents depend linearly 
on the order p. As mentioned in the introduction, the 
anomalous scaling (the deviation from linear dependence 
on p) of the exponents observed in many turbulent flows 
is a signature of deviations from scale invariance and of 
intermittency. The velocity and helicity scaling expo- 
nents in the direct cascade range of runs A and B are 
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FIG. 8; Probability density functions at different intervals in the direct cascade for velocity (solid) and helicity (dashed) 
increments in the direction perpendicular to the axis of rotation. Increments are normalized by their variance. The dotted 
curve represents a Gaussian distribution with the same variance. 



shown in Fig. [T] 

In run A, with Vl = 0.06, the effect of rotation is negli- 
gible and the velocity scaling exponents display the usual 
deviation from the Kolmogorov p/3 scaling. Deviations 
from a the straight line are often quantified in terms of 
the intermittency coefficient fx = 2(^3 — Ce^ which for this 
run is /i = 0.29 ± 0.06, in agreement with previous simu- 
lations and experiments of non-rotating turbulence. The 
third order exponent is ^3 = 0.99 ± 0.03, in good agree- 
ment with the value of 1 expected for isotropic and homo- 
geneous turbulence. The higher orders computed in this 
run are also consistent with results of non-rotating tur- 
bulence at very large Reynolds numbers (see e.g., [i^). 
All values of the scaling exponents up to order 8 are given 
in Table U 

On the other hand, for run B at low Rossby number, 
the velocity scaling exponents are (within error bars) con- 
sistent with a scale invariant (intermittency-free) linear 
relationship. For this run, /i = 0.0 ±0.1, compatible with 
a value of /i = which corresponds to a completely scale 
invariant flow. A similar result was reported in an exper- 
imental study (although associated in that case with 
a linear scaling (^p = p/2 and in a range of scales that 



may correspond to an inverse cascade of energy). The 
values of for run B contrast with results obtained for 
the scaling exponents in non-helical rotating turbulence 

3 




42| and from labora- 
, where a reduction 



from DNS for the direct cascade 
tory experiments for decaying flows 
of the intermittency was observed but anomalous scaling 
as a signature of intermittency was still present; as an 
example, in Ref. [i^ , for a non-helical flow at late times 
with Ro « 0.07 it was found that ^ = 0.24 ± 0.02. The 
possibility that the different behaviors reported in the 
experiments can be ascribed to helicity is tantalizing. 

In the case of the helicity exponents ^p, the highest 
orders could not be measured for the reasons discussed 
above: the higher order dependence on the fields of the 
helical structure functions, and the associated stronger 
fluctuations observed. However, results for run A are in 
agreement with previous studies of the direct cascade of 
helicity in isotropic and homogeneous turbulence, with 
the helicity being more intermittent than the velocity 
field (note however differences in the definitions used here 
and in the analysis in [l^l). This is illustrated by the 
smaller values of ^p when compared with (p in run A 
(see Fig. [7] and Table [l|. On the other hand, results in 




FIG. 9: (Color online) Slices of the energy density (top left), 
vorticity intensity (top right), z component of the velocity 
(bottom left), and helicity density (bottom right), in run B 
aX t K, 30. Note the imprint of small scales in the vorticity 
and helicity (right column). 



run B at low Rossby number seem to be of a different na- 
ture. The exponents up to p = 4 are within error bars 
consistent with a hnear (self-similar) scaling w 0.73p 
(i.e., with a slightly larger slope than for C,p in the same 
run, see Fig. [7]), but ^5 departs from such scahng (see also 
Table [J). This departure would indicate intermittency in 
the helicity, a property that will be confirmed in the next 
section studying the probability density functions of the 
increments. 
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FIG. 10: (Color online) Three dimensional rendering of the z 
component of the velocity in the entire domain in run B at 
t ~ 30 (above), and a zoom on a subregion (below) showing 
the z component of the velocity in a column-like structure 
(left) and its helicity density (right). 



V. PROBABILITY DENSITY FUNCTIONS 

The identification of multi-fractal (as opposed to scale 
independent) scaling in turbulent flows based on scaling 
exponents present several difficulties. It is well known 
that transients and finite size effects can give spurious 
multi- fractal scaling [i^, and that logarithmic or sub- 
leading corrections to the spectrum can also lead to the 
same result [i^l- In this section we consider probability 
density functions of the velocity and helicity increments. 
In a scale invariant fiow, the velocity increments are ex- 
pected to be Gaussian. On the other hand, in an inter- 
mittent flow, probability density functions are expected 
to have strong non-Gaussian tails. 

In Fig. [8] we show the probability density functions 
(PDFs) of the velocity and helicity increments in the 
direction perpendicular to the axis of rotation for four 
values of = 0.4, 0.2, 0.06, and 0.02. The PDFs are 



normalized by their variance, and a Gaussian with the 
same variance is shown as a reference. The PDFs of ve- 
locity increments for t\_ = 0.4 and 0.2 are close to Gaus- 
sian. Note that these increments correspond to scales in 
the anisotropic inertial range (see Figs. [3] and [4]) . PDFs 
close to Gaussian were also observed for Su for other in- 
crements in this range. However, for £± = 0.4 a "bump" 
(which also decays as a Gaussian) can be observed in the 
tails. This bump disappears as smaller increments are 
considered, and seems to be associated with the forcing 
(note that the mechanical forcing acts between £ « 0.78 
and 0.90). 

The bump is more evident in the PDF of helicity in- 
crements with £± = 0.4. For both values £± = 0.4 and 
0.2 (in the direct cascade inertial range) the PDFs of 
helicity increments are different than for velocity incre- 
ments: deviations from Gaussianity are evident, and the 
PDFs show strong tails. The non-Gaussian tails are the 
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signature of the presence of strong gradients, and of in- 
termittency in the spatial distribution of heUcity. These 
results confirm - independently of the scaling exponents 
analyzed in the preceding section - that the anisotropic 
direct cascade of energy towards smaller scales is close to 
Gaussian and scale independent, while the direct cascade 
of hclicity is intermittent. 

For increments in the range of scales where Sp{£j_} 
is of the same order as Sp{i\\) for p > 4 (see e.g.. Fig. 
S]) , both the PDFs of velocity and of hclicity increments 
show non-Gaussian tails (see Fig. [8] for i± — 0.06 and 
0.02). This further confirms that at very small scales, 
the strongest events (in the tails of the PDFs) tend to- 
ward a recovery of isotropy, which would lead in turn to 
a traditional (and intermittent) direct cascade of energy. 
However, simulations at larger resolution and at different 
Reynolds and Rossby numbers will be required to verify 
if this recovery of isotropy occurs always near the dissi- 
pation range, or if a second isotropic inertial range de- 
velops when enough scale separation is available at scales 
smaller than the Ozmidov scale when rotation can pre- 
sumably be neglected. 



VI. STRUCTURES 

The results discussed in Sect. IIVI and IVl point towards 
an anisotropic and scale-invariant energy distribution in 
an intermediate range of scales (smaller than the forc- 
ing scale), and a highly intermittent hclicity distribu- 
tion in the same range of scales. In helical rotating tur- 
bulence, energy cascades both towards large and small 
scales, while helicity cascades towards small scales domi- 
nating the direct cascade inertial range. In light of these 
facts, it is of interest to look at the structures that arise 
in the flow. 

In Figure [9] arc displayed horizontal slices (in the plane 
perpendicular to the rotation axis) of the energy density, 
the vorticity intensity, the z component of the velocity, 
and the helicity density in run B with Ro = 0.06, at 
t w 30. While the energy density seems organized in 
large-scale patches (substantially larger than the forcing 
scale which corresponds roughly to 1/7 of the box), the 
vorticity intensity and helicity density show small-scale 
structures. This is consistent with an inverse cascade 
of energy and a direct cascade dominated by the helic- 
ity. Note however that the energy density distribution in 
space is different from the distribution observed in two di- 
mensional turbulence, where an inverse cascade of energy 
also takes place; in other words, the inverse cascade here 
differs from the purely two-dimensional case, be it only 
because the conservation of helicity (in the ideal case) in- 
duces the flow to keep some trace of three-dimensionality 
(and isotropy) at small scale, as noted before. Moreover, 
smooth structures at some intermediate scale can be ob- 
served in the helicity and the vorticity (see e.g. the left 
side of the box); these regions are also correlated with 
similar regions in the z component of the velocity. 



When run B is started from a previously isotropic 
state, the flow first becomes anisotropic and then a self- 
organization process starts that leads to the formation 
of columns. Those columns can be identified when the 
energy density, helicity density, or the vorticity inten- 
sity are visualized. The columns have strong Uz (see e.g. 
Fig. nop . However, a few columns can be distinguished 
from the rest, in that they have a strong updraft veloc- 
ity and concentrate in their core positive helicity, with 
strong relative helicity (strong alignment between veloc- 
ity and vorticity). These columns are stable, and we were 
able to track these columns in the simulation for over ten 
turnover times. 

Far from these structures, the fiow displays a myriad 
of small scales, as illustrated by the small-scale filaments 
in the hclicity density. These filaments are also organized 
in columnar structures, but the thick columns with net 
helicity live for much longer times. As a result, in real 
space the inverse cascade of energy can be identified as 
the system evolves in time and these helical structures 
merge with columns with a lesser amount of relative he- 
licity, increasing the characteristic width of the column 
when its energy density is visualized, but keeping the 
thickness of the helical core in the column approximately 
constant. This is the result of the helicity injected in the 
flow cascading directly to smaller scales, which allows for 
a localized helical column, but prevents the formation 
of a thick distribution of helicity in a column fllling all 
space. Accompanying the direct cascade, strong fluctu- 
ations of helicity are observed in the turbulent columns 
with the characteristic size of the vortex filaments. 



VII. CONCLUSIONS 

The analysis of the structure functions and of the prob- 
ability density functions of velocity and helicity incre- 
ments from data stemming from direct numerical simu- 
lations of helical rotating turbulence at high resolution, 
showed that, at least for the strongest events in the small 
scales, isotropy is recovered at sufficiently small scales. 
This is observed both in the collapse of the parallel and 
perpendicular structure functions, as well as in the tran- 
sition from near-Gaussian statistics of the velocity incre- 
ments in the inertial range towards PDFs with strong 
tails in the same range of scales as observed in the struc- 
ture functions. More studies will be required to see if this 
transition takes place for all orders when the Reynolds 
number is large enough, leading to a sufficient scale sep- 
aration. 

Concerning scaling exponents, wc also confirmed that 
in the anisotropic direct cascade range the velocity incre- 
ments are (within error bars) scale invariant (i.e., non- 
intermittent) while helicity increments are intermittent. 
This is further confirmed by the probability density func- 
tions, which show strong non-Gaussian tails for the helic- 
ity, and are near Gaussian for velocity increments in the 
anisotropic range. The scaling exponents for the hclicity 
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are consistent (within error bars) with a scale invariant 
dependence ~ 0.73p up to p = 4, but ^5 departs from 
a straight line. More data will be required to confirm a 
possible bi-fractal or multi-fractal scaling for the helic- 
ity, although we would like to point out that a bi-fractal 
scaling would be consistent with the two types of heli- 
cal structures observed in visualizations of the flow (the 
large-scale laminar columns and the small-scale vortex 
filaments) . and reminiscent of the behavior or other sys- 
tems with two type of structures, e.g., the Burgers equa- 
tion, which develops smooth ramps connected by sharp 



shocks. 
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